Thermal photon spectra and HBT correlations at full RHIC energy 
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Using a model for the evolution of the fireball created in Au-Au collisions at full RHIC energy 
which reproduces the hadronic single particle spectra and two-particle correlations we calculate 
thermal photon emission from the hot partonic and hadronic matter. We present predictions for 
both the emitted photon transverse momentum spectrum and Hanbury-Brown Twiss (HBT) corre- 
lations. Surprisingly, we find that due to the strong transverse flow in the model photon emission 
from hadronic matter is sizeable in the momentum range expected to be dominated by thermal 
emission and that consequently the HBT radius parameters do not reflect exclusively properties 
of the partonic evolution phase. We compare to another model calculation where pre-equilibrium 
photon emission is dominant to illustrate the differences of both approaches. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 

In the hot and dense system created in a heavy-ion col- 
lision, the relevant momentum scales for scattering pro- 
cesses drop as a function of proper time as a large number 
of particles is created and the system evolves towards 
equilibrium. Assuming that an expanding thermalized 
system is created after some initial timescale To, the rele- 
vant momentum scales inside the thermalized matter are 
then set by the temperature, which in turn drops as the 
fireball volume expands with time. Thus, by selecting 
a certain window in momentum for scattering processes, 
one is preferentially sensitive to a certain proper time 
window of the evolution. If this window is chosen such 
that it is below the scales associated with non-equilibrium 
initial scattering processes, one is, for (p) ks 3Tj with Tj 
the initial temperature of the system, in principle sen- 
sitive to the initial hot and dense phase where hadronic 
matter presumably undergoes a phase transition to the 
quark-gluon plasma (QGP) as seen in lattice QCD sim- 
ulations (see e.g. yL|). 

Since electromagnetic probes (such as real photons and 
dileptons) can escape from the system at all times due 
to the relative smallness of the electromagnetic coupling 
constant, a measurement of hard real photon emission 
can be expected to provide information about the early 
phases of fireball evolution and ultimately a window to 
observe the QGP directly HQ. Several such investiga- 
tions have been carried out for the CERN SPS (see e.g. 
S ) I n a schematic way this was addressed also in |(j 
where data measured by the WA98 data Q could indeed 
be used to find constraints for the equilibration time. In 
addition, the measured photon spectrum allowed to dis- 
tinguish between a boost-invariant longitudinal expan- 
sion and a scenario involving initial compression and re- 
expansion, the latter being favoured by other observables 
as well [|. 

In a recent paper, an approach has been suggested aim- 
ing at determining the relevant scales of the spacetime 
evolution of a fireball created in an Au-Au collision at 
full RHIC energy ■ It is the aim of this paper to ex- 



plore what amount of thermal photon emission is pre- 
dicted from the resulting evolution and to what degree a 
precise measurement of the photon spectrum is able to 
further constrain the evolution dynamics of the system. 
In addition to photon transverse momentum spectra, we 
also discuss predictions for Hanbury-Brown Twiss (HBT) 
correlation measurements of hard photons and their ca- 
pability to provide additional information about the early 
stages of the system. 



II. THE MODEL 

In this section, we present a model framework for the 
fireball expansion which in essence is a parametrization 
inspired by a hydrodynamical evolution of the collision 
system. This model has been shown to give a good de- 
scription of both single particle spectra and two particle 
correlations at RHIC simultaneously for a breakup tem- 
perature well below the phase transition temperature. It 
is described in greater detail in 0, here we only repeat 
the essential facts: 

For the entropy density at a given proper time we make 
the ansatz 



i(T,r} s ,r) = NR(r,r) ■ H(t) s ,t) 



(1) 



with t the proper time measured in a frame co-moving 
with a given volume element, i] s = | ln(|^) the space- 
time rapidity and R{r, t), H(t] Si t) two functions describ- 
ing the shape of the distribution and N a normalization 
factor. We use Woods-Saxon distributions 



R(r, t) = 1/11 + cxp 
H(t] s ,t) = 1/ (l + cxp 



Rc(t) 



Tj s - H c (t) 



(2) 



for the shapes. Thus, the ingredients of the model are 
the skin thickness parameters d ws and ?7 WS and the para- 
metrizations of the expansion of the spatial extensions 
R c {t), H c (t) as a function of proper time. From the dis- 
tribution of entropy density, the thermodynamics can 



2 



be inferred via the EoS and particle emission is then 
calculated using the Cooper-Frye formula. For simplic- 
ity, we apply the model at the moment only to cen- 
tral and almost central collisions and assume that the 
flow is built up by a constant acceleration aj_, hence 
Rc{t) = R^ + ^t 2 with R° c an initial radial extension as 
found in overlap calculations In |9(, the model parame- 
ters have been adjusted such that the model gives a good 
description of the data. 

III. THE PHOTON EMISSION RATE 

The complete calculation of the photon emission rate 
from a partonic medium to order a * h as been a very 
involved task El El El El El El ill The relevant 
processes include qq annihilation, QCD compton scat- 
tering, bremsstrahlung and qq annihilation with subse- 
quent scattering. For the present calculation, we use the 
parametrization of the rate given in El] . 
Vector mesons play an important role for the emission 
of photons from a hot hadronic gas. The first calcu- 
lation of such processes has been performed in |l7j in 
the framework of an effective Lagrangian. It has been 
found that the dominant processes are pion annihilation, 
7r + 7r~ — > /?7, 'Compton scattering', n^p — > 7r ± 7 and p 
decay, p — > 7r + 7r~7. 

Several more refined approaches have been made since 
then (for an overview, see [3). In the following, we will 
use a parametrization of the rate from a hot hadronic gas 
taken from [lfij . 

A. The photon spectrum 

The spectrum of emitted photons can be found by folding 
the rate with the fireball evolution. In order to account 
for flow, the energy of a photon emitted with momentum 
k^ = (k t , kt, 0) has to be evaluated in the local rest frame 
of matter, using E = k^u^ with UiArj s ,r,T) the local 
flow profile. Following the results in |3] we assume n = 
const. ■ n s and n± = const. ■ r with n± the transverse 
rapidity in the following. The distribution of entropy 
density is manifest in the dependence of the temperature 
T = T(j] s ,t,t) on the spacetime position. In order to 
account for the breakup of the system once a temperature 
Tp is reached, a factor 8(T — Tp) has to be included into 
the folding integral. 

In order to be able to compare with SPS data and our 
previous model calculations, we present the differential 
emission spectrum into the midrapidity slice y = 0. 

IV. PHOTON SPECTRA 

The resulting photon spectra are shown in Fig.^ I n ( a ) ; 
we compare the calculation for RHIC with a model cal- 
culation for SPS done in the framework outlined above 



(essentially this reproduces the results obtained in Q for 
SPS in a similar, but slightly less sophisticated frame- 
work) and the data obtained for SPS 158 AGeV 10% 
central Pb-Pb collisions 0. 

As expected, the RHIC photon spectrum is character- 
ized by a larger absolute normalization (reflecting a larger 
four-volume of radiating matter) and a less steep falloff, 
reflecting a higher temperature of the emission region 
and/or a larger amount of transverse flow. Roughly, we 
find a factor two more photons at 2 GeV at RHIC than 
at SPS. 

In order to investigate this ambiguity between flow and 
temperature in more detail, we compare in (b) and (c) 
the amount of photon radiation coming from the QGP 
phase and the hadronic phase, respectively. Surprisingly, 
while at SPS photon emission from the hadronic phase is 
a relatively small contribution and the measured photon 
spectrum directly reflects the initial QGP phase of the 
evolution and hence the initial temperature, this is not 
so at RHIC where the QGP and hadronic contributions 
are found to be equally important. 

In order to identify the reasons for this surprising be- 
haviour, we artificially switch off the transverse flow and 
recalculate the photon spectrum. The result is shown 
in Fig. (d). This clearly demonstrates that it is the 
comparatively strong flow in the late hadronic evolution 
at RHIC (which is in this framework motivated by the 
strong falloff of the HBT correlations with transverse mo- 
mentum, see |2j) which leads to the observed importance 
of hadronic contribution. 

An additional factor is the different cooling in the 
hadronic phase: the system created at SPS conditions 
is characterized by a large net baryon density at midra- 
pidity. This in turn implies more production of heavy 
baryonic resonances at the phase transition at SPS as 
compared to RHIC. Decay processes of these resonances 
in turn lead to an overpopulation of pion phase space as 
compared to equilibrium, more so at SPS than at RHIC, 
which in turn causes more efficient cooling (see H and 
for details). In essence, there is a longer-lasting hadronic 
phase at RHIC which is hotter on average, leading to 
more photon emission and increasing the relative contri- 
bution from the hadronic phase. 

This is, as it stands, rather unfortunate since the slope of 
the measured photon spectrum is no longer a good mea- 
sure for the initial temperature at RHIC but rather for 
the flow in the late phases, which is comparatively well- 
known from the measured hadronic spectra and HBT cor- 
relations already. 

V. HBT CORRELATIONS FOR PHOTONS 

Using the folding of the rate with the fireball evolution 
as emission function S(x, K) (describing the amount of 
photons with momentum K 11 emitted at spacetime point 
x^) we calculate the HBT parameters as |27ll2^| 

R 2 eide = (f) R 2 out = ((x-/3jf) R loas = (z 2 ) (3) 





FIG. 1: Upper left: The direct photon spectrum measured at the CERN SPS as compared with the model calculations for 
SPS (dotted) and RHIC (solid). Upper right: The measured spectrum at SPS, shown with the model calculation for SPS (solid) 
and the contributions from QGP (thin) and hadron gas (dotted). Lower left: The photon spectrum in the model calculation 
for RHIC (solid) and the contributions from QGP (thin) and hadron gas (dotted). Lower right: As before, but without taking 
transverse flow into account when evaluating the photon emission from a given volume element. 



with = Xfj 



(x^) and 



</(*)>(*) 



/ d 4 xf{x)S(x 1 K) 
fd*xS(x,K) 



(4) 



Since photons escape without rescattering from all space- 
time points, we expect i? s ide to be smaller than in a mea- 
surement of hadronic correlations, reflecting the fact that 
photons from the whole volume enter the emission func- 
tion at all times (in contrast to hadron emission which 
takes place from the surface until final breakup), and we 
expect the difference of i? ut and i? S idc to be more pro- 
nounced, based on the observation that photon emission 
is not dominated by a small time interval at the final 
breakup of the fireball but shows a broader distribution 
(see FigEJ). 

In Figs|3 ^ and O (left panels) we show the three corre- 
lation radii and compare with the hadronic HBT in the 
measured region |9| (we do not make any calculations for 
low momentum photons since the emission rates based on 
perturbative expansions cannot be expected to be valid 



in this regime). 

These expectations are to some degree confirmed. i? s ide 
appears to be unexpectedly larger for photons than for 
pions, but in order to deduce the geometrical size, one 
has to take into account the effects of flow, and here the 
less steep falloff of the photon i? S id c as compared to the 
7r + measurement (reflecting the fact that a sizeable num- 
ber of photons is emitted early before transverse flow 
is built up) seems to extrapolate to a smaller value at 
k t = 0. There is a pronounced difference between i? ou t 
and i? S idc as expected before, and i?i on g appears to be 
approximately the same for photons and pions. This is 
not a surprise since the question of the transverse emis- 
sion coordinate and hence the difference between surface 
and volume emission is irrelevant here. 
In order to gain more insight into the underlying physics, 
we present the photonic correlation parameters as they 
would appear if one could separate hadronic and QGP 
emission experimentally (Figs|31 01 and [5J right panels). 
Consistently, the hadronic correlation radii show larger 
absolute normalization and steeper falloff with fc t , both 
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FIG. 2: The emitted multiplicity of photons (solid line) and neg. pions (dotted line) at midrapidity in the momentum regions 
displayed in Figs. [3] 2] and as a function of proper time r. Pion emission from the fireball grows roughly with the surface 
area until finally the hot matter decouples and the remaining pions are freed, therefore the emission is strongly peaked towards 
late times. In contrast, photon emission takes place from the whole volume at all times and is mainly a function of the average 
temperature at a given time, hence there is a large contribution from early times. 
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FIG. 3: Left panel: 7? s ide as a function of transverse pair momentum kt for photons (solid) and ir + (dotted). Right panel: 
Decomposition of the photonic -Rsidc (solid) into correlations among photons emitted from the QGP (dotted) and hadron gas 
(dashed). 



features caused by the flow which at the same time ex- 
pands the geometrical radius and tends to descrease the 
size of the observed region of homogeneity with k t . 

In contrast, the early emission of QGP photons shows a 
comparatively small-sized system characterized by little 
flow. 

Note that i? S ide for the total emission is to a good ap- 
proximation an average of the individual contributions 
weighted by their magnitude. However, this does not 
hold for i? ou t which includes a temporal component as 
well: The sum of emission durations of hadronic and 
QGP phase has to be larger than the individual emis- 
sion durations. 



VI. COMPARISON WITH OTHER WORKS 

Photon interferometr y ha s been investigated in a number 
of studies [U HIH 1H |24| . Here, we focus on a recent 
study |13 based on the parton cascade VNI/BMS 
and a subsequent boost-invarinat hydrodynamical evolu- 
tion. 

This particular model shows two pronounced differences 
to the present approach — is contains a pre-equilibrium 
contribution to the photon yield (not included here) and 
it is based on a boost-invariant expansion whereas our 
approach is characterized by accelerated longitudinal ex- 
pansion. 

In order to make a meaningful comparison, we have to 
correct for the fact that the correlation parameters in 
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k t [GeV] k t [GeV] 

FIG. 4: i?out as a function of transverse pair momentum k t for photons (solid) and tt + (dotted). Right panel: Decomposition 
of the photonic -Rout (solid) into correlations among photons emitted from the QGP (dotted) and hadron gas (dashed). 
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FIG. 5: Ri ong as a function of transverse pair momentum kt for photons (solid) and ir + (dotted). Right panel: Decomposition 
of the photonic Ri on g (solid) into correlations among photons emitted from the QGP (dotted) and hadron gas (dashed). 



[25j are calculated for central collisions whereas in our 
approach we have chosen to calculate for the 30% most 
central collisions in order to be able to compare with the 
ir + HBT results of the model (which describe in turn the 
measurements). This amounts roughly to a 15% increase 
of the radii when we go from the values shown in Fias l3l4l 
and [5] to central collisions. 

In |25|. i? S idc at kt = 1 GeV is found to be 4.6 fm, at k t = 
2 GeV 3.8 fm. In our model we find for central collisions 
similar (though slightly smaller) values of 4.3 fm and 3.5 
fm respectively. Thus, i? s id for photons mainly reflects 
the initial spatial extension of the emission source and 
shows little signs of the late source expansion. 
The values for R out given in [25| are 4.7 fm for k t = 1 GeV 
and 3.7 fm for k t — 2 GeV, in contrast to 5.11 fm and 
3.9 fm in the present model. Presumably this is a conse- 
quence of the different distribution of emission in time: 
Whereas in [2{| the emission at kt = 2 GeV is almost 
completely dominated by pre-equilibrium contributions 



with extremely short duration (leading to i? ut ~ -Rsidc) 
this is different in our model, causing -R ou t > -Rsidc- 
^?iong again shows only a moderate difference between 
the two models: In 25], at kt = 2 GeV a value of 1.6 fm 
is found as compared to 1.95 fm in the present model. 
In [25j, this is a consequence of the dominance of pre- 
equilibrium photons in this momentum region which are 
mainly emitted at times r < 0.3 fm/c, giving little room 
for longitudial expansion, the observed value of -Riong is 
then mainly dictated by the uncertainty principle. In 
contrast, the emission of thermal photons in the present 
model starts only at 0.6 fm/c when the system has under- 
gone some expansion. Moreover, early longitudinal flow 
is stronger for a boost-invariant expansion, leading to a 
further reduction of .Riong in 25] . Given the fact that we 
do not include a pre-equilibrium component of photon 
emission, our value of 2 fm probably overestimates the 
data somewhat. 

The effect of a pre-equilibrium contribution on the HBT 
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correlations depends on the relative magnitude of the 
contribution at given transverse momentum. Currently, 
there is still considerable uncertainty about the absolute 
magnitude of a pre-equilibrium contribution caused by 
questions like the role of LPM suppression. According to 
|29j this may reduce the pre-equilibrium yield by up to 
a factor 3.5 at 2 GeV transverse momentum. Presum- 
ably, a measurement of i?i on g would be most capable of 
answering the question if a pre-equilibrium contribution 
constitutes a sizeable fraction of the measured photon 
spectrum below 3 GeV or not. 

VII. SUMMARY 

We have presented predictions for the spectrum of ther- 
mal photons at full RHIC energy and for two photon 
correlation radii. This was done in a model which is able 
to describe both the measured hadronic single-particle 
spectra and two-particle correlations. 
We find that the photon spectrum between transverse 
momenta 2 GeV < kt < 4 GeV is not dominated by 
emission from the hot QGP phase as for SPS conditions 
[f|. Instead, thermal emission from the hadronic phase 
turnes out to be a significant contribution. This is mainly 
due to the strone transverse flow at RHIC (which in turn 
is required for a good description of the hadronic HBT 



correlations). 

Photon HBT correlations are expected to show a smaller 
system characterized by less flow than seen in hadronic 
HBT measurements. This is a consequence of surface ver- 
sus volume emission processes. However, in the present 
approach HBT correlations are not exclusively domi- 
nated by the QGP phase but reflect an average of the 
QGP and hadronic evolution phases. Given the im- 
portance of thermal photon radiation from the hadronic 
phase, this is not surprising. 

If this scenario turns out to be realized at RHIC then an 
extraction of QGP properties from photon spectra and 
correlations may prove difficult. Possibly the cleanest 
window would then be the high momentum tail domi- 
nated by pre-equilibrium physics where the number of 
emitted photons may still help to understand how fre- 
quently scattering processes occur. 
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